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ABSTRACT 



Context. A study of late-type low-mass eclipsing binaries provides us with important information about the most common stars in the 
Universe. 

Aims. We obtain the first light curves and perform period analyses of two neglected eclipsing binaries GK Boo and AE For to reveal 
their basic physical properties. 

Methods. We performed both a period analysis of the times of the minima and a BVR light curve analysis. Many new times of minima 
for both the systems were derived and collected from the data obtained by automatic and robotic telescopes. This allowed us to study 
the long-term period changes in these systems for the first time. From the light curve analysis, we derived the first rough estimates of 
the physical properties of these systems. 

Results. We find that the analyzed systems are somewhat similar to each other. Both contain low-mass components of similar types, 
both are close to the Sun, both have short orbital period, and both contain another low-mass companions on longer orbits of a few 
years. In the case of GK Boo, both components are probably of K3 spectral type, while the distant companion is probably a late M 
star. The light curve of GK Boo is asymmetric, which probably causes the shift in the secondary minima in the O-C diagram. System 
AE For comprises two K7 stars, and the third body is a possible brown dwarf with a minimal mass of only about 47 M Jup . 
Conclusions. We succeed in completing period and light curve analyses of both systems, although a more detailed spectroscopic 
analysis is needed to confirm the physical parameters of the components to a higher accuracy. 

Key words, binaries: eclipsing - stars: fundamental parameters - stars: individual: GK Boo, AE For 
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1. Introduction 

Low-mass sta rs are the most common stars in our Galaxy (e.g. 
lKrounal l2002>. However, owing to their low luminosity, only 
these close to the Sun have been studied in detail and many of 
them have never been analyzed. Hence, we focused on two rather 
neglected low-mass eclipsing binary systems: GK Boo and AE 
For. Their light curves as well as their period modulation had 
never been studied. Some stu dies indicate that most late-type 
stars are single (e.g. Lada 2006), but the number of papers study- 
ing the multiplicity of the late-type systems is still rather limited. 
Therefore, the incidence of multiples in late-type stars remains 
unexplored. 

The study of eclipsing binaries provide us with important 
information about the physical properties of both of their com- 
ponents - their radii, masses, and evolutionary status. However, 
when considering only with the light curve, several assumptions 
have to be made. For the analysis presented in this paper we 
also used the photometric data obtained by automatic and robotic 
telescopes (such as ASAS, Pi of the sky, and SWASP). Thanks 
to these huge databases of observations, the long-term evolution 
of these systems can be studied for the first time. 



2. GK Boo 

2.1. Introduction 

The system GK Boo (= BD+37 2556, V max = 10.86 mag) is 
an Algol-type eclipsing binary with an orbital period of about 
0.48 day. It is also a primary component of a visual double 
designated WPS J14384+3632 in t he Washington Double Star 
Catalog (WDSQ, Mason et al. 2001). The secondary component 
of this double star is about 14" distant, and is probably gravita- 
tionally bound to GK Boo itself. It is about 0.4 mag fainter, but 
since its discovery in 1933 there has been no detectable mutual 
motion of the pair, hence the orbital period is of about thousands 
of years (rough estimation from the Kepler's law). 

The star is too faint, thus was not observed by Hipparcos 
satellite, and its distance is therefore rather uncertain. 
iKharchenkol 1200 1 ) introduced the parallax 30.29 mas, which is 
however only an estimate. Its spectral type is also un known, but 
the B — V index derived from the Tycho catalogue (H0g et al., 
2000), B — V — 0.89 mag indicates a spectral type of about Kl. 
On th e other hand, the 2MASS infrared photometry (ICutri et al.L 
l200l g ives J—H = 0.527 mag (therefore a spectral type of K3). 
Finally, lAmmons et al.l (2006) introduced a temperature corre- 
sponding to a spectral type of about K2-3. Al l the se rough spe c- 
tral estimates were taken from lPopperl d!980l) and |Coi (2000). 



* This paper uses observations made at the South African 
Astronomical Observatory (SAAO) 



http://ad.usno.navy.mil/wds/ 
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2.2. Light curve 



The star was observed by the SuperWASP dPollacco et all f2006) 
project and its complete light curve (hereafter LC) is available. 
However, we did not use these data for the LC analysis be- 
cause these were not measured in any standard photometric fil- 
ter. These data were only used to derive the minima times (see 
below). We observed the target at the Ondfejov observatory in 
the Czech Republic with the 65-cm telescope equipped with the 
CCD camera. For the light curve analysis, only the data from 
two nights in May 2011 were used (see the electronic data ta- 
bles). The remaining observations were used for the minima 
time derivation and to analyze the period changes in the system 
(see below section l2~3l l. The observations were obtained in stan- 
dard B , V, and R filters according to the specification of lBesselll 
(fl990l) . 

At firs t, the complete LC w as analyzed using the program 
PHOEBE dPrsa & Zwitterll2005l). which is based on th e Wilson- 
Devinney algorithm (WD. I Wilson & Devinnevl[l97ll) . The de- 
rived quantities are as follows: the secondary temperature T2, 
the inclination i, the luminosities L,, the gravity darkening co- 
efficients gi, the albedo coefficients A,, and the synchronicity 
parameters Fj. The limb darkening was approximated using 
a linear law, and the values of x, were interpo lated from the 
van Hamme's tables, given in Ivan Hammel(ll993h . 

At the beginning of the fitting process, we fixed the tempera- 
ture of the prima ry compone nt at T\ = 4700 K (corresponding to 
spectral type K3, Cox 2000). In the absence of spectroscopy, the 
mass ratio was derived via a so-called "q-search method". This 
means that we tried different values of mass ratio in the range 1 .5 
- 0.5 in steps of 0. 1 and tried to find the best LC fit according to 
the lowest value of rms. Finally, we found that the best-fit solu- 
tion was reached with the value q = Mz/M\ = 0.9, which agrees 
with both eclipses having almost equal depths. For a given mass 
ratio, the semi-major axis was fixed to an appropriate value for 
the primary mass to be equal to a typical mass o f a particular 
spectr al type (e.g. Popperi ll980l lHarmanedll98"8l or lAndersenl 
Il99ll) . With this approach, we were able to estimate the masses, 
in addition to the radii of both components in absolute units. 

However, during the LC fitting process we found that the 
LC of GK Boo is asymmetric. In particular, the part of the LC 
near the secondary minimum is distorted in all BVR filters. The 
brightness just after the ascent from the secondary minimum 
(near the phase 0.6) is higher than the brightness just before 
the descent (phase 0.4). The difference is about 0.022 mag in 
B, 0.018 mag in V, and 0.017 mag in R filter, respectively. 

With the PHOEBE code, we tried to fix the values of A, 
and gi to their appropriate values of 0.5 and 0.32, respectively. 
However, after then we also allowed these parameters to be fit- 
ted, because the fit is tighter (rms). However, probably owing 
to the asymmetry of the LC these quantities converged to the 
rather improbable values given in Table [1] and the shape of the 
observed LC could not be fitted properly. For the asymmetry of 
the curve, we also tried to introduce a star spot on either of the 
components. However, no acceptable solution with spot(s) was 
found to describe the shape of the light curve more accurately in 
the PHOEBE program. The parameters of the LC fit are given 
in TableQ] but these cannot sufficiently describe the shape of the 
LC. 

We therefore tried a different code, called ROCHE, devel- 
oped by Theo Pribulla (Pri bullaL 120041) . which is also based 
on the WD code but has for instance also some other comput- 
ing methods and different controlling of the calculation process. 
With this program, we used two star spots and similar input pa- 




Fig. 1. Light curves in BVR filters for GK Boo, the solid lines 
represent the final fit. The residuals after the fit are plotted below. 
The curves are shifted along y-axis for reasons of clarity. 



rameters as described above. At the beginning, the values of A,- 
and gi were fixed to the appropriate values of 0.5 and 0.32, re- 
spectively. However, to achieve a tighter fit both A, and gi values 
were also varied across the range from to 1 in steps of 0.05 for 
both components. The synchronicity parameters converged to 
much more reliable values. The value of mass ratio was fixed 
to q = 1.0 and then also fitted as a free parameter. This was 
possible because there is a clear distorti on of the LC outside the 
minima (see e.g. lTerrell & W ilson 2005). For the fitting process, 
the two different limb darkening laws were also tried, namely a 
linear and logarithmic. The latter one provides a much tighter fit 
to our data. All of the resulting LC parameters are also given in 
Table Q] (together with parameters of two cooler spots located on 
the primary component - longitude, latitude, radius and temper- 
ature factor). As one can see, the two solutions clearly differ even 
outside their respective error bars for some of the parameters. 

The individual errors in the parameters were not taken from 
the WD code, but derived in the following way. We computed a 
range of solutions for GK Boo, which were then used for its error 
estimation. All solutions with x 1 value close to the minimal one 
(5% from our final solution) were taken and the resultant values 
of parameters were used to compute the differences between the 
parameters. The errors in the individual parameters were then 
computed as a maximum difference and their individual WD er- 
rors, given by maxicii — a m j„) + 5ai + 6a m i„. 

This solution obtained with the ROCHE program provides a 
much closer fit to the observed data and is the fit plotted in Fig. 
Q] The value of the eccentricity was fixed at (for a discussion 
about possible eccentricity see below). Our resultant parameters 
indicate that both the components are still located on the main 
sequence, (as required because the age of the Universe does not 
allow low-mass stars to have evolved from the main sequence). 
If we follow the assumption of a K3V primary, then the sec- 
ondary is also of K3V spectral type. These are consistent with 
the photometric indices presented above, as well as with the in- 
dividu al masses and radii for these types of stars (e.g. lHarm anec 
1988). An undetectable value of the third light was also resulted 
derived by this analysis. The presence of photospheric spots on 
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Table 1. Light curve parameters of GK Boo. 



Parameter 


Value 




PHOEBE 


ROCHE 


r, [K] 


4700* 


r 2 [K] 


4540 ± 50 


4615 ± 63 


q(=M 2 /M 1 ) 


0.9±0.1 


0.95 ±0.12 


e 


0* 




i [deg] 


89.83 ± 0.57 


89.28 ± 0.37 


gi 


0.00 ± 0.04 


0.35 ± 0.05 


g2 


0.00 ± 0.03 


0.35 ± 0.05 


A, 


0.00 ± 0.08 


0.80 ± 0.05 


A 2 


1.00 ± 0.08 


0.80 ± 0.05 


F, 


1.892 ±0.107 


1.131 ± 0.096 


Fi 


1.866 ±0.116 


1.295 ±0.108 


U (B) [%] 


54.8 ± 1.9 


52.4 ± 1.1 


L 2 (B) [%] 


45.2 ± 1.8 


47.6 ± 1.0 


L, (V) [%] 


53.5 ± 1.5 


51.6 ± 1.2 


L 2 (V) [%] 


46.5 ± 1.3 


48.4 ± 1.1 


U (R) [%] 


52.1 ± 1.4 


51.1 ± 1.1 


L 2 (R) [%] 


47.9 ± 1.3 


48.9 ± 1.0 




Spots: 




li [deg] 


— 


287.2 ± 7.9 


bi [deg] 


- 


60.5 ± 3.2 


ri [deg] 




37.9 ± 2.0 


ki 




0.75 ± 0.04 


1 2 [deg] 




63.3 ±7.4 


b 2 [deg] 




47.4 ± 12.8 


r 2 [deg] 




28.7 ±4.1 


k 2 




0.76 ± 0.04 




Derived quantities: 




*1 [Ro] 


0.83 ±0.18 


0.89 ±0.15 


#2 [Rg] 


0.86 ±0.18 


0.86 ±0.14 


M, [Mo] 


0.73 ± 0.06 


0.73 ± 0.06 


M, [Mo] 


0.66 ± 0.06 


0.70 ± 0.06 


Note: * - fixed. 



both components of such a late spectral type star is also foresee- 
able. 

2.3. Period analysis 

To monitor the detailed long-term evolution of the sys- 
tem or its short-period modulation, we collected all avail- 
able pub lished minima obs ervations. Photometry from the 
SWASP dPollacco et all |2006l) ASAS dPoimanskil 12002). and 
PiOfTheSky dBurd et al.L 120051) projects were used to derive 
many new minima times for GK Boo. All of these data are 
given in Table [TUl which is available in el ectronic form only. 
The method of Kw ee & van Woerdenl (1956) was used. Some of 
the data were of poor quality, but most were accurate enough to 
perform a detailed period analysis of the system. The range of 
these data is about 12 years. 

We used these data to analyze the period modulation and 
found some interesting results. Applying the hypothesis of a 
third body in the system (the so-called Li ght- Time Effect, here- 
after LITE, described e.g. bv lIrwinlfT959h . we found a weak pe- 
riod modulation with a period of about four years. The final fit 
to the data together with the theoretical curve is shown in Fig. 
|2] As one can see, there is also some long-term period evolution 
of the orbital period (the blue dashed line), which was described 
as a quadratic term in ephemerides. It can be understood as a 
slow period decrease caused by the mass loss from the system 
or mass flow between the components (or even momentum loss, 
magnetic breaking, etc.). Another explanation is that this is only 



Table 2. Final parameters of the long orbit for GK Boo. 



Parameter 


Value 


HJD 


2454305.4570 ± 0.0006 


P [day] 


0.47777174 ± 0.00000022 


Pi [day] 


1472.7 ± 170.0 


Pi [yr] 


4.032 ± 0.450 


A [day] 


0.0126 ± 0.0012 


To 


2454263.3 ± 1108.3 


<jl> 3 [deg] 


56.54 ± 15.0 


^3 


0.084 ± 0.267 


£[-10- 10 ] 


-1.071 ± 0.206 


f(M 3 ) [M Q ] 


0.000633 ± 0.000002 


M X min [Mq] 


0.115 ±0.001 


M x( ,o [M Q ] 


0.134 ±0.002 


M x -io [M ] 


0.242 ± 0.005 


a l2 sin i [AU] 


0.217 ± 0.108 


a 3 [mas] 


88.7 ±9.8 



part of the long-term period modulation, although we have only 
limited data coverage. 

A more interesting finding is that of a period of about 4 years. 
Applying the LITE hypothesis, we obtained a final set of param- 
eters given in Table [2j namely the period of the third body 773, 
the semi-amplitude of the effect A, the time of periastron pas- 
sage Tq, the argument of periastron 0)3, and the eccentricity e^. 
Despite the low amplitude (about only 1.8 minutes) of the LITE, 
most of the observed minima times are of higher precision and 
the modulation is clearly visible. Table [2] also provides the mass 
function of the third body /(M3), which helps us to estimate its 
predicted mass. 

Having no information about the inclination between the or- 
bits of the eclipsing pair and the hypothetical third body, we 
plotted Fig. [3] where a plot mass versus inclination is shown. 
Assuming the coplanar orbits (i.e. z'3 = 90° — > M3 = M3 m! „), 
the resulted minimum mass of the third body is only about 
0.116 M Q , which places this body at the lower end of stellar 
masses, hence we can rule out the hypothesis of a brown dwarf 
or even an exoplanet. Despite of there being no upper limit to 
this mass (it goes to infinity with 23 — > 0°), we can estimate a 
lower limit to the mass. Taking into account that no third light 
is detected in the LC solution, e.g. £3/(^4 + L2) < 0.01 and 
assuming a main-sequence star, we can estimate its mass to be 
lower than 0.22 M , which is shown in Fig. [3] as a gray area. 
Further observations are still needed to confirm this hypothesis 
with higher conclusiveness. 

If we assume the parallax of GK Boo as given by 
Khar chenkol(l200ll) . n — 30.29 mas, we are also able to compute 
the angular distance of a hypothetical body to be about 89 mas. 
This separation of components is well above the limit for mod- 
ern stellar interferometers. However, there is a problem with the 
brightness of the third component, which was found to be about 
more than five magnitudes fainter than the eclipsing pair itself. 
With the brightness of about 1 1 mag for the system, this makes a 
detection impossible. The magnitude difference of the third body 
with respect to the close pair also clarify why no third light was 
detected in the LC solution. 

Another interesting result was a detection of displaced sec- 
ondaries. This can be clearly seen in more precise data points 
(SWASP and our new observations). That secondary minima oc- 
cur at a different phase of <p2 + 0.5 from the primary usually 
indicates that the system is on an eccentric orbit. GK Boo is a 
well-detached system, so the eccentric orbit cannot be ruled-out 
easily. Therefore, we assumed an apsidal motion hypothesis for 
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Table 4. Methods of minima fitting for GK Boo. 
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Fig. 2. Periodic modulation of period GK Boo. Blue dashed line 
represents quadratic ephemeris, while red solid line stands for 
the LITE hypothesis. Residuals are plotted in the bottom plot. 
The larger the symbol, the higher the weight (higher the preci- 
sion). 
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Fig. 3. GK Boo: Mass of the third body based on from the LITE 
hypothesis with respect to the inclination between the orbits. 

Table 3. Apsidal motion parameters for GK Boo. 



Parameter 


Value 


e 


0.0944 + 0.0068 


co [deg] 


268.9 ± 2.5 


6j [deg/cycle] 


0.00026 ± 0.00001 


U [yr] 


1790 ± 50 



our da ta set of minima times. We followe d a procedure described 
by e.g . iGimenez & Garcfa-Pelavol dl983h or lGimenez & Basterol 
( 1995) and obtained a set of apsidal motion parameters. The plot 
of residuals (after subtraction of the LITE fit) with the apsidal 
motion fit is shown in Fig. [4] It is obviously very slow because 
the position of secondaries versus primaries changes only very 
slowly. The resultant values of apsidal motion parameters are 
given in Table [3] 

However, we have to rule out this hypothesis because it lead 
to unacceptable results. With some information about the physi- 
cal parameters of both components, we can use the apsidal mo- 
tion parameters to estimate the internal s tructure con stant. The 
theoretical log k.2,theor value taken from IClaretl d2004l) should 
range from -1.35 to -1.65. However, the mean value of log k2 of 
both components that can be derived from our solution is very 
different, even when k2 < 0, which is unacceptable. Thus, the 
system is very probably on a circular orbit. 

We may ask why the secondary minima deviate from the 0.5 
phase. We published a finding that the displaced secondary min- 
ima can also be present i n contact binaries where no eccentric 
orbit is possible dZaschel 1201 ll) . so one cannot perform an ap- 
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0.00119 
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Fig. 4. O - C diagram of GK Boo with the apsidal motion hy- 
pothesis, black color is for primary minima, while blue one for 
secondary. 



sidal motion analysis based only on the minima times of a par- 
ticular system. Some studies found that the secondary minimum 
is displaced because of the distortion of the LC, thus any stan- 
dard routine for deriving the time of minimum (e.g. Kwee-van 
Woerden, bisector chord method or polynomial fitting) cannot be 
used properly because these consider symmetric minima only. 
When using these methods to determine minima where both as- 
cending and descending branches have different slopes, we re- 
cover only a "false eccentricity". 

One can also ask about a significance of the fits presented 
in Figs. [2] and [4] For this comparison, we summarized different 
approaches in Table[4] In addition to the rms values, we also pro- 
vide th e values of BIC (Bayesian Information Criterion, see e.g. 
lLiddldl2007h . which show the significance of the fit. According 
to this method, the smaller the rms value, the tighter the fit. To 
conclude, our final fit provides the smallest rms, but its signifi- 
cance is low and still highly speculative. This is also caused by 
the poor data coverage, and large scatter in the minima and their 
low accuracy. Determinations of new more precise minima are 
therefore needed to confirm or exclude this hypothesis. 



3. AE For 

3.1. Introduction 

The Algol-type system AE For (= HIP 14568, V max = 
10.22 mag) is also a poorly studied binary. Its published spectral 
types range from K4 to M0, with the m ost probable one being 
K7V as derived by dTorres et all 120061) . The system was pre- 
sented as a wide double with the star HD 19632 based on the ir 
similar parallaxes and proper motions (see lPoveda et al.l [T994). 

Neither the light curve nor the radial velocity curve of 
AE For have been studied. The star was observed by the 
Hipparcos satellite and a few times also for the minima obser- 
vations. It was also continuously monitored with automatic pho- 
tometric systems such as PiOfTheSky and AS AS. However, the 
quality of these data do not allow us to use them for a LC anal- 
ysis. The distance to the system was derived from the Hipparcos 
data to be d — 3 1 .5 pc. 
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Table 5. Light curve parameters of AE For. 
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Fig. 5. Light curve in BVR filters for AE For, the solid line rep- 
resents the final fit. The curves are shifted along the y-axis for 
greater clarity. 



3.2. Light curve 

We observed the star from the South African Astronomical 
Observatory (SAAO) in 2010, using the classical one-channel 
photoelectric photometer mounted on the 50-cm telescope. All 
measurements w ere carefully reduced to the Cousins E-region 
standard system dMenzies et al.L Q~989) and corrected for differ- 
ential extinction. 

Thanks to its orbital period close to one day, its complete 
light curve was observed once in standard BVR filters, with some 
overlapping points (about 170 data points in each filter were ob- 
tained). Unfortunately, the quality of the data acquired for sev- 
eral nights was not very good, hence the scatter in the curve is 
affected by these conditions. Two secondary and one primary 
minima were observed (see below). 

We analyzed our data using the same computational proce- 
dure as for GK Boo. The primary temperature was fixed to the 
appropriate value of 4100 K (sp K7V), the eccentricity was fixed 
to 0, the values of gravity darkening coefficients were fixed at 
0.32, and the albedo coefficients to 0.5 (as recommended for 
stars with convective envelopes), while the limb d arkening co- 
efficie nts were interpolated from values given in Ivan Ha mme 
( 19931). The computational approach was different for the mass 
ratio q, which was fixed to q = 1 .0 because of the weak outside- 
eclipse ellipsoidal variations and its detached configuration. In 
addition, the synchronicity parameters Fj were set to values of 
1.0 for both components. The program ROCHE was used and 
the resulting LC parameters are given in Table [5] while the final 
solution is presented in Fig. [5] 

One can see that the secondary temperature T2 is close to the 
value of T[, indicating that the components are similar. Thus, the 
estimated spectral types of both stars are probably K7V + K7 V. 
Both components are still located on the main sequence and their 
properties are in agre ement with the ty pical values of K7V stars 
(as presented by e.g. lHarmaned ll988). The third light was also 
not detected here in any filter. In contrast to GK Boo, the LC of 
AE For seems to be symmetric. 



Parameter 


Value 


r, [K] 


4100* 


r 2 [K] 


4065 ± 48 


q(=M 2 /Mi) 


1.0* 


e 


0* 


i [deg] 


86.51 +0.31 


Si = g2 


0.32* 


Ai = A 2 


0.50* 


F { =F 2 


1.000* 


L, (B) [%] 


63.2 ± 1.3 


L 2 (B) [%] 


36.8 + 1.0 


Li (V) [%] 


63.1 ± 1.2 


L 2 (V) [%] 


36.9 ± 1.0 


Li (R) [%] 


62.6 ± 1.4 


L 2 (R) [%] 


37.4 ± 1.0 


Derived quantities: 


Ri [R Q ] 


0.66 ±0.10 


Ri [Ro] 


0.52 ± 0.08 


Af, PVLJ 


0.50 ± 0.05 


Mi [Mo] 


0.50 ± 0.05 


Note: * - fixed. 
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Fig. 6. O - C diagram of AE For. Up: With linear ephemeris. 
Bottom: With respect to the phase. The data points are fitted with 
the curve representing the third body hypothesis (see the text for 
details). 



3.3. Period analysis 

Similarly to GK Boo, we tried to perform the period analysis of 
all available minima. The collection of minima i s much smaller, 
but th anks to the first observation by Hipparcos dPerrvman et all 
fl997l) these cover a longer time span than for GK Boo. Several 
new minima were derived based on our new observations from 
SAAO as well as those from the ASAS and PiOfTheSky surveys. 

The same hypothesis as for GK Boo was applied to the data 
points here. All of the minima times used for the analysis are 
summarized in Table QT| which is available in electronic form 
only. As one can see from Fig. [6] there is a clear variation in the 
minima times. We used the same third-body hypothesis (LITE) 
as for GK Boo, deriving a final fit to the data given by the pa- 
rameters in Table[6] The LITE hypothesis resulted in a rather ec- 
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Table 6. Final parameters of the long orbit for AE For. 



Parameter 


Value 


HJDq 


2452605.97070 ± 0.00035 


P [day] 


0.91820943 ±0.00000012 


P3 [day] 


2524.6 ± 149.6 


Pi [yr] 


6.912 ± 0.409 


A [day] 


0.00083 ± 0.00032 


To 


2453548.8 ±413.1 


a>3 [deg] 


146.2 ± 57.8 


£3 


0.601 ±0.414 


f(M 3 ) [M ] 


0.000098 ± 0.000001 


M X min [Mo] 


0.047 ± 0.001 


M 3 , 60 [Mq] 


0.055 ± 0.001 


^3,30 [Mo] 


0.098 ± 0.003 


a\2 sin i [AU] 


0.167 ±0.064 


a 3 [mas] 


117.2 ± 8.3 



Table 7. Methods of minima fitting for AE For. 



Method of minima fitting 


rms 


BIC 


Linear ephemeris: 


0.000255 


24.4 


LITE and linear ephemeris: 


0.000163 


44.8 



centric orbit, although the result is affected by a relatively large 
error, hence maybe the value should be lower. Only additional 
observations would help us confirm or refute this hypothesis, re- 
fine the period, and possibly detect some long-term evolution of 
the period similar to that in GK Boo, because the first obser- 
vation from Hipparcos deviates significantly from the fit. With 
the same procedure as for GK Boo, we computed the signifi- 
cance of the fits according to the BIC criterion (see Table [7]). 
As one can see, the fit is still very poor and highly speculative. 
However, using only the linear ephemeris, there remains a clear 
quasi-sinusoidal variation, which needs some physical explana- 
tion. 

From the LITE parameters, we were able to calculate the 
minimal mass of the third body (i.e. coplanar orbits), which 
we found to be only about 47 Mj up , which is even lower than 
the limit of stellar masses. Therefore, if the orbits were copla- 
nar (which only would be our assumption, because the pro- 
cess of tidal coplanarization is very slow), the third body would 
very probably be a brown dwarf (exoplanets have masses about 
one half of this value). With such a body, we reach minimal 
masses that can be detected by this method, because the ampli- 
tude of LITE is comparable to the typical precision of individual 
minima-time measurements. Whatever applies to the possible in- 
terferometric detection of GK Boo companion also applies here, 
because its luminosity is too low. 

4. Discussion and conclusions 

We have derived preliminary light-curve solutions and period 
analyses of the poorly studied Algol-type eclipsing binaries 
GK Boo and AE For, which we have found to have several in- 
teresting and similar features. Since both of them are low-mass 
stars of very similar types (K3+K3 for GK Boo, and K7+K7 for 
AE For), both of them have short orbital periods. Moreover, both 
are relatively close to the Sun and also appear to contain third 
bodies in their systems, which cause a periodic modulation of 
the orbital periods of both systems. Assuming a coplanar orbit, 
for AE For this third body appears to be a brown dwarf, which 
makes this system even more interesting. However, more photo- 



metric and spectroscopic observations are needed to confirm or 
refute this hypothesis. 

The system GK Boo has an asymmetric light curve, which 
is the probably accounts for the shift in the secondary minimum 
in phase with the primary one. The apsidal motion hypothesis 
cannot explain this discrepancy. 

In general, if the third body hypothesis as proposed based 
on the period analysis is found to be the correct one, here we 
have considered quite curious examples of hierarchical quadru- 
ple systems of low masses. As far as we know, there are only 
a few similar multiple late-type systems for which one of the 
components is an eclipsing binary (e.g. BB Scl or MR Del). 
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Table 8. Heliocentric minima of GK Boo used for the analysis. Table 9. Minima of GK Boo, cont. 



HJD - 2400000 


Error 


Type 


Filter Reference 


HJD - 2400000 


Error 


Type 


Filter 


Reference 


51260.8547 


0.0007 


sec 


IBVS 5060 


54163.55942 


0.00065 


prim 


- 


SWASP 


51273.9890 




prim 


A.Paschke - Rotse 


54165.70909 


0.00037 


sec 


- 


SWASP 


51311.7377 


0.0004 


prim 


IBVS 5060 


54166.66462 


0.00022 


sec 


- 


SWASP 


53128.46588 


0.00024 


sec 


SWASP 


54167.62053 


0.00034 


sec 


- 


SWASP 


53128.70498 


0.00041 


prim 


SWASP 


54170.72638 


0.00023 


prim 


- 


SWASP 


53129.66141 


0.00056 


prim 


SWASP 


54171.68178 


0.00043 


prim 


- 


SWASP 


53130.61697 


0.00047 


prim 


SWASP 


54189.59755 


0.00032 


sec 


- 


SWASP 


53132.52777 


0.00021 


prim 


SWASP 


54190.55322 


0.00019 


sec 


- 


SWASP 


53135.39435 


0.00010 


prim 


SWASP 


54191.50702 


0.00053 


sec 


- 


SWASP 


53137.54360 


0.00024 


sec 


SWASP 


54194.61418 


0.00029 


prim 


- 


SWASP 


53138.49891 


0.00028 


sec 


SWASP 


54195.57050 


0.00029 


prim 


- 


SWASP 


53139.45491 


0.00021 


sec 


SWASP 


54202.49743 


0.00046 


sec 


- 


SWASP 


53141.60560 


0.00027 


prim 


SWASP 


54206.55909 


0.00039 


prim 


- 


SWASP 


53152.59424 


0.00025 


prim 


SWASP 


54208.46922 


0.00027 


prim 


- 


SWASP 


53153.54982 


0.00020 


prim 


SWASP 


54208.70823 


0.00143 


sec 


- 


SWASP 


53154.50553 


0.00017 


prim 


SWASP 


54210.61908 


0.00030 


sec 


- 


SWASP 


53155.46071 


0.00015 


prim 


SWASP 


54212.53056 


0.00050 


sec 


- 


SWASP 


53156.41642 


0.00014 


prim 


SWASP 


54213.48698 


0.00049 


sec 


- 


SWASP 


53158.56492 


0.00048 


sec 


SWASP 


54214.44267 


0.00029 


sec 


- 


SWASP 


53159.52074 


0.00022 


sec 


SWASP 


54214.68147 


0.00039 


prim 


- 


SWASP 


53160.47640 


0.00021 


sec 


SWASP 


54215.63620 


0.00057 


prim 


- 


SWASP 


53161.43169 


0.00026 


sec 


SWASP 


54216.59233 


0.00019 


prim 


- 


SWASP 


53163.58313 


0.00031 


prim 


SWASP 


54217.54805 


0.00038 


prim 


- 


SWASP 


53164.53862 


0.00031 


prim 


SWASP 


54218.50359 


0.00021 


prim 


- 


SWASP 


53165.49392 


0.00015 


prim 


SWASP 


54219.45928 


0.00030 


prim 


- 


SWASP 


53166.44965 


0.00020 


prim 


SWASP 


54220.41494 


0.00069 


prim 


BVR 


PS 


53169.55361 


0.00034 


sec 


SWASP 


54220.41498 


0.00029 


prim 


- 


SWASP 


53170.50939 


0.00024 


sec 


SWASP 


54222.56423 


0.00066 


sec 


BVR 


PS 


53171.46489 


0.00034 


sec 


SWASP 


54222.56444 


0.00028 


sec 


- 


SWASP 


53172.42077 


0.00031 


sec 


SWASP 


54223.51949 


0.00169 


sec 


BVR 


PS 


53175.52707 


0.00024 


prim 


SWASP 


54223.51994 


0.00108 


sec 


- 


SWASP 


53176.48289 


0.00047 


prim 


SWASP 


54224.47565 


0.00035 


sec 


- 


SWASP 


53180.54262 


0.00022 


sec 


SWASP 


54225.43110 


0.00022 


sec 


- 


SWASP 


53181.49835 


0.00017 


sec 


SWASP 


54225.67010 


0.00026 


prim 


- 


SWASP 


53182.45376 


0.00018 


sec 


SWASP 


54226.38671 


0.00034 


sec 


- 


SWASP 


53183.41009 


0.00076 


sec 


SWASP 


54226.62544 


0.00036 


prim 


- 


SWASP 


53191.53267 


0.00080 


sec 


SWASP 


54227.58111 


0.00026 


prim 


- 


SWASP 


53192.48678 


0.00031 


sec 


SWASP 


54228.53676 


0.00031 


prim 


- 


SWASP 


53193.44381 


0.00026 


sec 


SWASP 


54230.44800 


0.00019 


prim 


- 


SWASP 


53194.39850 


0.00045 


sec 


SWASP 


54231.40365 


0.00019 


prim 


- 


SWASP 


53197.50383 


0.00080 


prim 


SWASP 


54231.64168 


0.00048 


sec 


- 


SWASP 


53198.46046 


0.00023 


prim 


SWASP 


54232.59790 


0.00032 


sec 


- 


SWASP 


53199.41623 


0.00026 


prim 


SWASP 


54233.55304 


0.00025 


sec 


- 


SWASP 


53203.47632 


0.00029 


sec 


SWASP 


54234.50899 


0.00034 


sec 


- 


SWASP 


53209.44918 


0.00047 


prim 


SWASP 


54235.46550 


0.00054 


sec 


- 


SWASP 


53215.42110 


0.00040 


sec 


SWASP 


54236.42047 


0.00028 


sec 


- 


SWASP 


53220.43588 


0.00083 


prim 


SWASP 


54236.65842 


0.00085 


prim 


- 


SWASP 


53221.39317 


0.00039 


prim 


SWASP 


54237.37517 


0.00137 


sec 


BVR 


PS 


53226.40796 


0.00058 


sec 


SWASP 


54239.52553 


0.00091 


prim 


BVR 


PS 


53231.42630 


0.00045 


prim 


SWASP 


54249.55879 


0.00027 


prim 


- 


SWASP 


53232.38275 


0.00032 


prim 


SWASP 


54250.51471 


0.00031 


prim 


- 


SWASP 


53237.39732 


0.00033 


sec 


SWASP 


54251.47023 


0.00057 


prim 


- 


SWASP 


53243.37030 


0.00082 


prim 


SWASP 


54252.42644 


0.00028 


prim 


- 


SWASP 


53248.38731 


0.00078 


sec 


SWASP 


54254.57647 


0.00059 


sec 


- 


SWASP 


53259.85637 


0.00078 


sec 


SWASP 


54256.48657 


0.00172 


sec 


BVR 


PS 


53827.68574 


0.00085 


prim 


SWASP 


54256.48655 


0.00039 


sec 


- 


SWASP 


53829.59740 


0.00088 


prim 


SWASP 


54257.44222 


0.00025 


sec 


- 


SWASP 


53830.55201 


0.00069 


prim 


SWASP 


54261.50367 


0.00031 


prim 


- 


SWASP 


53831.50729 


0.00051 


prim 


SWASP 


54262.45860 


0.00020 


prim 


- 


SWASP 


53832.46270 


0.00015 


prim 


SWASP 


54263.41487 


0.00037 


prim 


- 


SWASP 


53832.70138 


0.00024 


sec 


SWASP 


54265.56380 


0.00079 


sec 


- 


SWASP 


53833.65681 


0.00047 


sec 


SWASP 


54266.51892 


0.00042 


sec 


- 


SWASP 


53837.47803 


0.00029 


sec 


SWASP 


54267.47491 


0.00020 


sec 


- 


SWASP 


53851.57305 


0.00105 


prim 


SWASP 


54268.43084 


0.00017 


sec 


- 


SWASP 


53852.52690 


0.00049 


prim 


SWASP 


54271.53663 


0.00070 


prim 


- 


SWASP 


53853.48213 


0.00027 


prim 


SWASP 


54272.49225 


0.00033 


prim 


- 


SWASP 


53854.43799 


0.00032 


prim 


SWASP 


54273.44808 


0.00026 


prim 


- 


SWASP 


53855.39312 


0.00090 


prim 


SWASP 


54276.55314 


0.00095 


sec 


- 


SWASP 


53855.63171 


0.00037 


sec 


SWASP 


54277.50813 


0.00065 


sec 


- 


SWASP 


53856.58695 


0.00051 


sec 


SWASP 


54278.46353 


0.00031 


sec 


- 


SWASP 


53887.40789 


0.00024 


prim 


SWASP 


54279.41942 


0.00021 


sec 


- 


SWASP 


53901.50285 


0.00080 


sec 


SWASP 


54305.45758 


0.00157 


prim 


BVR 


PS 


54140.62657 


0.00022 


prim 


SWASP 


54307.37045 


0.00141 


prim 


BVR 


PS 


54149.70458 


0.00021 


prim 


SWASP 


54568.70649 


0.00072 


prim 


" 


Piofthesky 


541 50 65998 


00078 




SWASP 


54568.94557 


0.00097 


sec 




Piofthesky 


54153.76487 


0.00119 


sec 


SWASP 


54925.36514 


0.0001 


sec 


R 


OEJV 107 


54154.71989 


0.00026 


sec 


SWASP 


54937.3112 


0.0020 


sec 


Ir 


IBVS 5918 


54155.67581 


0.00061 


sec 


SWASP 


54937.5462 


0.0005 


prim 


Ir 


IBVS 5918 


54156.63161 


0.00053 


sec 


SWASP 


54947.34260 


0.0014 


sec 


R 


OEJV 107 


54157.58667 


0.00065 


sec 


SWASP 


54958.8085 


0.0007 


sec 


V 


IBVS 5894 


54159.73799 


0.00019 


prim 


SWASP 


54959.52489 


0.0001 


prim 


R 


OEJV 107 


54160.69347 


0.00018 


prim 


SWASP 


55354.40380 


0.00148 


sec 


BVR 


PS 


54161.64928 


0.00022 


prim 


SWASP 


55364.4367 


0.0001 


sec 


BVRI 


IBVS 5965 


54162.60325 


0.00019 


prim 


SWASP 


55385.45891 


0.00131 


sec 


BVR 


PS 
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Table 10. Minima of GK Boo, cont. 



HJD - 2400000 


Error 


Type 


Filter 


Ref. 


55386.41381 


0.00104 


sec 


BVR 


PS 


55391.43075 


0.00097 


prim 


BVR 


PS 


55392.38658 


0.00113 


prim 


BVR 


PS 


55590.66335 


0.00036 


prim 


- 


RU 


55599.50146 


0.00018 


sec 


- 


RU 


55616.46275 


0.00005 


prim 


B 


PZ 


55616.46266 


0.00012 


prim 


R 


RU 


55619.56796 


0.00010 


sec 


R 


RU 


55634.61898 


0.00016 


prim 


R 


RU 


55640.58947 


0.00019 


sec 


R 


RU 


55644.41169 


0.00005 


sec 


B 


PZ 


55650.62303 


0.00014 


sec 


R 


RU 


55651.34031 


0.00113 


prim 


B 


PZ 


55662.56706 


0.00013 


sec 


R 


RU 


55671.40664 


0.00007 


prim 


VR 


PZ 


55685.50021 


0.00012 


sec 


R 


RU 


55687.41150 


0.0001 1 


sec 


BVR 


PZ 


55692.42853 


0.00004 


prim 


BVR 


PZ 


55700.55071 


0.00004 


prim 


BVR 


PZ 


55707.47805 


0.00008 


sec 


BVR 


PZ 



Table 11. Heliocentric minima of AE For used for the analysis. 



HJD - 2400000 


Error 


Type 


Filter 


Reference 


48500.6581 


0.001 


prim 


Hp 


Hipparcos 


51180.9089 




prim 


R 


VSOLJ 37 


51180.9090 


0.0002 


prim 


R 


VSOLJ 47 


51191.0099 




prim 


R 


VSOLJ 37 


51191.0100 


0.0002 


prim 


R 


VSOLJ 47 


51191.9279 




prim 


V 


VSOLJ 37 


51191.9280 


0.0001 


prim 


V 


VSOLJ 47 


51196.9770 


0.0001 


sec 


R 


VSOLJ 47 


51196.9774 




sec 


R 


VSOLJ 37 


5 1 504 1 1 886 

-) 1 JUt - . 1 lOOU 






I 


VSOT T 37 


51504.1190 


0.0002 


prim 


I 


VSOLJ 47 


52065.14178 


0.00121 


prim 


V 


ASAS 


52065.60355 


0.0025 


sec 


V 


ASAS 


52235.0140 




prim 


I 


VSOLJ 39 


52240.9825 




sec 


I 


VSOLJ 39 


52258.8860 


0.0002 


prim 


R 


VSOLJ 39 


52605.9710 




prim 


I 


VSOLJ 40 


52818.07823 


0.00077 


prim 


V 


ASAS 


52818.53514 


0.00165 


sec 


V 


ASAS 


52901.1754 




sec 


I 


VSOLJ 42 


52929.1801 




prim 


I 


VSOLJ 42 


52936.0674 




sec 


I 


VSOLJ 42 


52957.1865 




sec 


V 


VSOLJ 42 


52970.0410 




sec 


I 


VSOLJ 42 


52976.0097 




prim 


I 


VSOLJ 42 


52987.0267 




prim 


V 


VSOLJ 42 


52987.9466 




prim 




VSOLJ 42 


53300.1379 




prim 


V 


VSOLJ 43 


53335.0300 




prim 


V 


VSOLJ 43 


53340.0796 




sec 


V 


VSOLJ 43 


53705.0677 


0.0002 


prim 


V 


VSOLJ 44 


53728.02219 


0.00093 


prim 


V 


ASAS 


53728.48316 


0.00194 


sec 


V 


ASAS 


54039.75501 


0.00027 


sec 




PiOfTheSky 


54040.67301 


0.00045 


sec 




PiOfTheSky 


54047.1000 


0.0001 


sec 


V 


VSOLJ 45 


54052.60840 


0.00149 


sec 




PiOfTheSky 


54448.81542 


0.00039 


prim 




PiOfTheSky 


54726.11607 


0.00187 


prim 


V 


ASAS 


54726.57421 


0.00131 


sec 


V 


ASAS 


54862.9297 


0.0002 


prim 


V 


VSOLJ 50 


55558.0139 




prim 


V 


VSOLJ 51 


55570.41065 


0.00072 


sec 


BVR 


PZ - SAAO 


55571.32812 


0.00105 


sec 


BVR 


PZ - SAAO 


55576.37847 


0.00065 


prim 


BVR 


PZ - SAAO 
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